Electronic states on a twin boundary of a <i-wave superconductor 
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We show that an induced s-wave harmonic in the superconducting gap of an orthorhombic d x 2_ y 2 
superconductor strongly affects the excitation spectrum near a twinning plane. In particular, it 
yields bound states of zero energy with areal density proportional to the relative weight of the 
s-wave component. An unusual scattering process responsible for the thermal conductivity across 
the twin boundary at low temperatures is also identified. 
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The controversy concerning s-wave versus <i-wave pair- 
ing in high-T c superconductors is being gradually re- 
solved in favor of the d x 2_ y 2 symmetry. This is possible 
because of precise measurements of spontaneously gener- 
ated half-integral flux quanta on bicrystal and tricrystal 
films jjj together with the corner SQUID experiments 0, 
which all suggest that the superconducting gap changes 
sign on the Fermi surface between the a- and 6-directions 
with zeros or nodes in between. The existence of corre- 
sponding low-energy excitations in the superconducting 
state of cuprates have also been demonstrated by mea- 
surements of the magnetic penetration depth ||] , the nu- 
clear spin relaxation rate Q, and by angular resolved 
photoemission data ||. 

However, not all experimental data fit into a simple 
picture of a pure d x 2_ y 2 order parameter. Sun et al. 
H have found a non-vanishing tunneling current along 
the c-axis between YBa2Ca306+ K (YBCO) and the con- 
ventional superconductor Pb, which cannot exist for a 
tetragonal d-wave superconductor. It was argued |7]]|] 
that such an experiment can be understood if an admix- 
ture of an s-wave component due to the orthorhombicity 
of YBCO crystals |9j is taken into account. The mixed 
wave d+s gap then has nodes on the Fermi surface shifted 
from the principal diagonal directions. Recent high reso- 
lution measurements of the a6-plane thermal conductiv- 
ity have failed to see this effect imposing an upper limit 
of 10% on the weight of the s-wave harmonic [10|. 

An important challenge for current research is to rec- 
oncile various conflicting results on symmetry of the su- 
perconducting order parameter in YBCO by direct de- 
termination of the s-wave component in the gap. In this 
work we study the structure of the quasiparticle spec- 
trum close to a twinning plane in a predominantly d-wave 
orthorhombic superconductor and show that it changes 
significantly in the presence of an s-wave component of 
the gap. In particular, there appear bound states of zero 
energy with areal density proportional to the magnitude 
of the s-wave component. Bound states cause a modi- 
fication of the local density of states which can be ob- 
served by scanning tunneling microscopy. Other new ef- 



fects studied below include an unusual low temperature 
mechanism of heat transport across the twin boundaries. 
This is a kind of Andreev transmission process because 
it involves particle-hole conversion. 

The tetragonal symmetry of the Cu02 planes in YBCO 
is spoiled by the presence of CuO chains. The two possi- 
ble orientations of the chains lead to the formation of so- 
called twin domains separated by twin boundaries (TBs), 
which are (110) or (110) planes of the original tetrago- 
nal lattice (Fig|l|). The structural anisotropy gives rise to 
an anisotropy in the pairing interaction and, hence, to an 
admixture (real combination) of s- and d-waves in the su- 
perconducting gap with the relative phase or n in differ- 
ent twins ||. Since we assume the dominant component 
of the gap to be of the d-wave symmetry, a natural energy 
requirement is that ipd stays almost constant across the 
boundary, whereas ip s has different signs on both sides, 
thus exhibiting a soliton-like structure. The symmetry 
of the combined order parameter (d ± s) is, in this case, 
odd with respect to the reflection in the twinning plane. 
Note, that the opposite assumption of an even symmetry 
state (s ± d) is in conflict with the afo-tunneling exper- 
iments on heavily twinned YBCO samples Sigrist 
and co-workers suggested that the s-wave component 
could change sign between two twin domains by avoiding 
zero at the twin boundary and forming a time-reversal 
breaking state d + se %x . Existence of the time-reversal 
breaking states on twin boundaries in YBCO is still an 
open question, we therefore leave discussion of the elec- 
tronic spectrum for such a twin boundary to the end. 

The soliton-like behavior of the induced s-wave compo- 
nent plays a key role in the formation of bound states on 
twin boundaries. Examples of fermions trapped on one- 
dimensional inhomogeneities of the off-diagonal potential 
have been known for a long time pj}| . To get some in- 
sight into this we first consider a very simplified model of 
the boundary between two twins. Namely, we neglect all 
orthorhombic effects except the induced inhomogeneous 
s-wave component and assume cylindrical Fermi surface. 
The Bogoliubov-de Gennes (BdG) equations in the qua- 
siclassical approximation p2[, i.e., retaining terms of the 
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lowest order in (kp^o) , are written for a slowly varying 



function of space coordinates (u, v) 
the usual BdG wave function, as 



' r ip, tp being 
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where v F = kp/m (Ti = 1). For the pair potential close 
to the twin boundary (x — plane) we choose 



A(k, r) = A d (k) + A s tanh(a;/^o) 



(2) 



which gives correct orthorhombic gaps A^ ± A s at x — > 
±oo. At E — 0, system (|l|) can be transformed to a pair 
of independent first-order differential equations, which 
are easily solved giving 
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where k' F = {-k Fx ,k Fy ), A d (k' F ) = -A d (k F ), £ d = 
vp x /\Ad(kp)\, £ s — v Fx /A s . These solutions correspond 
to physical states whenever normalization condition can 
be satisfied. Since the asymptotic behavior of i[> does not 
depend on £qi we immediately conclude that the bound 
state for a given k F exists if A(kp, oo) and A(kp, — oo) 
have different signs, e.g., for k y c < k Fy < k y A in Fig. |^. 
This requirement can also be written as |Ad(kp)| < A s . 
For a model <i-wave gap Ad(k) = A (k^ — kf)/k F or, 
in our coordinate system, A<z(k) = 2A<jk x k y /k F , the 
midgap excitations exist in the vicinity of (±k F ,0) for 
\k y \ < (A s /2Ao)fc_F and close to (0,±fc F ) with the same 
condition on k x , though in the latter case, as we shall 
see below, these states are easily destroyed by various 
perturbations. 

The local density of states near the twin boundary has 
a (5-like peak at E = corresponding to the bound states. 
The associated areal density 
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is proportional to the relative weight of the s-wave com- 
ponent. Furthermore, it is possible to show with the help 
of the Atiyah-Singer theorem that the zero-energy branch 
found for the gap function (||) exists for a wide class of 
gaps between two twins |13[| . 

To describe more realistically anisotropy of electron 
properties in YBCO we now consider a second model, 
which is a tight-binding model of a single CuC>2 plane. 
The model includes the orthorhombicity through unequal 
effective hopping matrix elements t\ and ti along the a 
and b axes. The two different hopping amplitudes allow 
us to construct a twinning plane on a square lattice as 
shown in Fig. [l]. An extra local potential (which we de- 
note by Uq and assume positive) associated with atoms 



on the twin boundary is allowed by the symmetry and 
will be shown to have important consequences for the 
twin boundary properties. Thus, the one-electron Hamil- 
tonian of the normal phase is 

h = - m%kW + u o c l c i° ■ ( 5 ) 

ija ieTBer 

The Schrodinger equation for the discrete function ip.^ is 



(6) 



The dispersions of plane waves in the two twins are 
e(k) = -2(ti <2 cos k a + t 2 .i cos k b ) - fx 

= —At cos —= cos —= =F Ate sin —= sin —= — fi , (7) 
v2 v2 v2 v2 

where ti 2 = t(l ± e), and the lattice constant is taken to 
be unity. 

Quasiparticles scatter on the twin boundary preserv- 
ing their energy and the parallel component of their mo- 
mentum. To study the scattering process we first find 
from Eq. (^) normal components of the Fermi momenta 
for incoming and outgoing electrons with a given parallel 
component k y (see Fig. ^|). In the right twin they are 
k± x — ±fc + q, while in the left twin k± x = ±fc — q, 
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k — ^cos(q/v / 2) 
V2 Atcos(k y /V2) 
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The normal components of the Fermi velocities for all 
these states have, however, the same absolute value v Fx = 
2\/2tsin(A:/v / 2)cos(fc y /V2)/cos(q/v / 2). (Note, that ve- 
locities are defined with dimension of energy.) 

The wave function describing the scattering of an elec- 
tron approaching the twin boundary from the right is 



■ r, +re !k +' r - (3?>0), ip?=de lK - Ti (x<0). (9) 



Two equations for the unknown amplitudes r and d are 
(i) continuity of the wave function ip^ — ip^ and (ii) 
Eq. (|^) for x — 0, which we write, after some algebra and 
Fourier transformation over y, as 



iq/V2 



/V 2 = Up cos(q/V2) 
2t cos(k y /V2) 



(10) 



Remarkably, if Uq — 0, quasiparticles move across the 
twinning plane without scattering for arbitrary degree of 
orthorhombicity e. The potential barrier on the bound- 
ary yields a nonzero probability of the reflection with 
amplitude r = ce/(i — a), a = Uq/ 1 \[2v Fx . 

Below T c we again will not attempt to solve the prob- 
lem self-consistently but rather use a "guess" order pa- 
rameter to study qualitative features of the quasiparticle 
spectrum. The discrete BdG equations are 
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where ip stands now for a vector with quasielectron and 
quasihole components. For orthorhombic symmetry and 
nearest-neighbor interactions we define Aa± a = Ai/2, 
A it ± b = -A 2 /2 in one twin and A ii±a = A 2 /2, A ii±b = 
-Ai/2 in the other, Ai )2 = A(l ± S) (see Fig. 0). The 
superconducting gap in the two twins is given by Fourier 
transform of the off-diagonal element in (11): 

A(k) = A(cos k a — cos ± A<5(cos k a + cos kb) ■ (12) 

A nonzero parameter 5 describes effect of crystal or- 
thorhombicity on superconducting properties and deter- 
mines a relative weight of the s-wave component. 

In our model of the twin boundary the first term of 
Eq. (p"2|), i.e., the e?-wave component, is chosen to be 
constant, whereas the second, extended s-wave part has 
a discontinues change of sign at the twinning plane. Al- 
though this assumption is not correct and A s varies on 
a characteristic scale £ s 3> a, the qualitative picture re- 
mains the same since the discussed results depend only on 
the asymptotes of the order parameter. In the following 
analysis we will use the symmetry of the problem coming 
from relations ts-i^j = Uj anc i A&i&j = —Ay, where a is 
reflection in the twinning plane. The BdG equations are 
invariant under the combined transformation U = tsct, 
T3 being the Pauli matrix which acts in the particle-hole 
space. All solutions of Eqs. ( |ll|) are, hence, classified by 
their parity with respect to U. 

In the bulk of each twin, Eqs. ( pr|) are diagonalized by 
a transformation to plane waves. The wave function of a 
bound state is their linear combination 



ik-L — K+X 



fi± = JA 2 ± - E 2 , k± = n±/v Fx 



(13) 



for x > 0, while on the left side ipf for even and odd 
eigenf unctions are obtained from ipf — ±r 3 ?/;? i . The 
energy of the bound state and the parameter R are de- 
termined from Eq. (|ll| ) for i 6 TB. In the leading or- 
der in A, that is, neglecting corrections 0(A 2 /ef) to 
the energy, the characteristic equation coincides with 
Eq. (|l(|) for the two-component function Jl3|) . Solving 
it for Uq = we find that all bound states have zero en- 
ergy independent of normal state anisotropy. They exist 
if the gaps of incident and outgoing quasiparticles satisfy 



A+A_ > , 



(14) 



e.g., for k y c < k y < k y A and for k y > k y T> and k y < k y B 
in Fig. p. Taking into account the symmetry of the twin 
boundary it is easy to see that this condition coincides 
with the previously obtained A(kp, oo)A(kp, — oo) < 0. 



It is, however, opposite to the condition on the existence 
of surface bound states in a e?-wave superconductor [jlJJ . 

The boundary potential Uq splits the degeneracy be- 
tween even and odd states. For quasiparticles moving 
nearly perpendicular to the twinning plane, the dimen- 
sionless parameter a = Uo/V^vfx ~ Uo/ef is small. 
Bound states exist in the same region as in unperturbed 
case ( |l4| ) and their energies acquire dispersion 
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even, odd 



2aA + A_ 



A, 



A J 



(15) 



which yields a finite width of zero-energy peak in the 
local density of states (^). The effect of the boundary 
potential is destructive for bound states with wave vec- 
tors k± nearly parallel to the twin boundary, i.e., for 
bound states with k y > k y D and k y < k y -Q in Fig. |^. 
A weak barrier Uq/ef ~ A s /Ao is sufficient to destroy 
completely these bound states. 

As a further check of the robustness of the midgap 
bound states, we have relaxed the condition of constant 
nearest-neighbors pairing amplitudes in the Hamiltonian 
and changed them near the twin boundary as shown in 
Fig. |l]. This change models a space dependent gap on the 
discrete lattice and again does not affect energies of the 
bound states irrespective of the difference between A's 
and A"s. 

The transmission of excitations by the twin bound- 
aries, which is important for the heat conduction process 
at low temperatures, also has unusual features. In par- 
ticular, for T < A s the heat is carried to the twin bound- 
ary by low energy excitations with k « kA in Fig. |^. For 
these quasiparticles |A_| < E < \ A + \ w A s . Let us con- 
sider such an electron-like state with k « k^ approaching 
the boundary from the right. It cannot be reflected back 



into the twin as an electron with k 



nor transmit- 



ted across the boundary as an electron with k « kf, 
since the gap magnitude for the outgoing quasiparticles 
is |A_|_| > E. However, this electron can transform into a 
hole having k « k^ and the same gap A_ but with the 
opposite group velocity (Andreev reflection |^5[). In the 
case of zero barrier (Uq — 0) the Andreev reflection is 
perfect (there is only an exponentially decaying wave in 
the left twin). The transmission of heat across the twin 
boundary must thus come from excitations with E > A s 
and there will be a temperature dependent contribution 
to the thermal conductivity from the twinning planes of 
the form ktb <x exp(— Ag/ksT). 

An extra potential Uq on the twin boundary, however, 
opens up an additional scattering channel in which an in- 
cident electron-like quasiparticle transforms into an out- 
going hole-like excitation with k ss k^ on the other side 
of the boundary (Fig. |||). It is this process of Andreev 
transmission, which will be responsible for the thermal 
conductivity across the twin boundaries in the low tem- 
perature limit. Interestingly, the boundary barrier Uq, 
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which causes reflection of the quasiparticles from the twin 
boundary in the normal state, allows transmission of the 
low-energy excitations in the superconducting state. 

Finally, we comment on properties of bound quasi- 
particles for the time-reversal breaking superconducting 
state of the twin boundary fij] . These effects can be most 
easily understood on the basis of Eq. ( fl|). We assume that 
the gap function has, in addition to (0), a small imagi- 
nary component of the s-wave harmonic iA s f(x), f{x) 
being an even function. The right-hand side of Eq. (jl|) 
contains in this case an extra term — A s T2/(x). First- 
order energy corrections found by averaging this pertur- 
bation with respect to the unperturbed states (§) are 
nonzero and have opposite signs for ip + and ?/>_. As a 
result, the bound state peak (|J) in the local density of 
states shifts to finite energies, making possible to distin- 
guish the two superconducting states on twinning planes 
in tunneling experiments. Condensation of bound states 
with a nonzero momentum produces also a local current 
flow parallel to the twin plane, which was derived in |?]] 
in the framework of the Ginzburg-Landau theory. 

In conclusion, the formation of electron bound states 
on twin boundaries is a general property of orthorhombic 
predominantly <i-wave superconductors. Observation of 
a bound state peak in the local density of states on twin- 
ning planes in YBCO by means of scanning tunneling 
microscopy would be a direct confirmation of the pres- 
ence of s-wave component in Cu02 planes. Also, a novel 
process of Andreev transmission is responsible for ther- 
mal conductivity across the twin boundaries in the low 
temperature limit. 
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